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Abstract 

This paper presents a new type of surface models constructed on the basis of Finsler 
geometry. A Finsler metric is defined on the surface by using an underlying vector 
field, which represents an in-plane tilt order. According to the orientation of the vec- 
tor field, the Finsler length becomes dependent on both position and direction on the 
surface, and for this reason the parameters such as the surface tension and bending 
rigidity become anisotropic. These parameters make the surface non-spherical. To 
confirm that the model is well-defined, we perform Monte Carlo simulations under 
several isotropic conditions such as those given by random vector fields. The results 
are comparable to those of previous simulations of the conventional model, which 
is isotropic. It is also confirmed that a tubular phase appears when the vector field 
is constant. 

Key words: Surface model, Anisotropic membranes, Phase transition, Finsler 
geometry 

PACS: 64.60.-i, 68.60.-p, 87.16.D- 



1 Introduction 

Membranes are understood as a mapping from a two-dimensional surface M to 
R 3 pp. If M is of sphere topology, the image X(M) corresponds to a spherical 
membrane. The shape of X(M) is governed by a surface tension energy and 
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a bending energy, which are described by a surface model of Helfrich and 
Polyakov (HP) [2]|3] . As a general framework, the Landau- Ginzburg (LG) free 
energy is assumed [4.5.6J. The LG free energy includes higher order terms of 
the field X. Owing to such mathematical backgrounds, the surface shape and 
its phase structure have been extensively studied [7|8f9|10|llj . 




An anisotropic phase is predicted in a surface model with a bending rigid- 
ity which is anisotropic in the internal direction of the surface [12713] . The 
existence of such anisotropic surface was numerically confirmed [H]. These 
anisotropic surface models have also been studied by non-perturbative renor- 
malization group formalization [15]. A tubular surface can also be seen in a 
surface model with elastic skeleton, where a one-dimensional bending energy 
is assumed [To] . 



However, the origin of the anisotropy in membranes still remains to be clarified. 
Indeed, it is possible to ask why the bending rigidity becomes isotropic (or 
anisotropic). It is well-known that there exist anisotropic membranes without 
skeletons [T7] . 



Therefore, it is interesting to study a surface model on the basis of Finsler 
geometry [18]. In Finsler geometry, an infinitesimal length unit is considered 
to be dependent on the directions, and it gives a more general framework 
than the Riemannian geometry. Because of this potential ability in Finsler 
geometry, several studies have been conducted on its applications to physics 



In this paper, a Finsler geometric (FG) surface model for membranes is stud- 
ied and Monte Carlo (MC) simulation data are presented. This model is con- 
structed by extending a discrete surface model of Helfrich and Polyakov such 
that the metric function is replaced by a Finsler metric. The assumed Finsler 
metric is defined by using a vector filed v on M. Since the Euclidean length 
of v depends on the direction, this vector field v is considered to be an ori- 
gin of surface anisotropy. We should note that the isotropy is restored if v 
is given locally at random. In this case, the FG model should have the same 
phase structure as the conventional surface model of Helfrich and Polyakov. 
We firstly check this under several conditions and confirm that the FG sur- 
face model is well-defined. Nextly, it is demonstrated that a tubular surface is 
obtained by assuming that v is constant. 
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2 Finsler geometric surface model 



2. 1 Elements of Finsler geometry 



In this subsection, we briefly summarize the elements of Finsler geometry 
Let M be a two-dimensional manifold, and let C be a curve on M such that 
C 3 t t- >■ x(t) G M. We call M as a Finsler space if there exists a Finsler 
funtion L on M such that the Finsler length s of the curve C is given by 

t 

L{x,y)dt, (1) 

to 

where L(x,y) is a homogeneous function of degree 1 with respect to y. The 
symbols x = (xi,x 2 ) and y = (2/1,2/2) — (dxi/dt,dx 2 /dt) in L are a point on 
C and a tangential vector at x with the direction along which t increases, 
respectively. Thus, we have 

L(x,ky) = kL(x,y) (2) 

for any positive k. This equation implies that the Finsler length s of the curve 
C is independent of the parameter t. For this reason, the Finsler length s 
depends only on the ratio y 2 /yi, because y = (2/1,2/2) — 2/1(1,2/2/2/1) can be 
replaced by (1, 2/2/2/1) i n Eq. fl2}. We should note that the definition of Eq. ([T]) 
can also be written as 

^ = ^y)- ( 3 ) 

An example of the Finsler function L(x(t),y(t)) is given by using a vector field 
v such that 

Hx{t),y{t)) = <lTtf/\v\, (4) 



where |v| = ^JJ2i{dxi/ds) 2 . We call y^iUi an d |v| as the Euclidean lengths. 
Note that the reparametrization dxi/dt— (dxi/ds)(ds/dt) allows us to write 



Y^iVi = yY^i(dxi/ds) 2 ds/dt. Thus, we have ds/dt = yJ2iVi /l v l- This leads 
to the expression s in Eq. (OQ). The vector v along C in Eq. (j4j) can also be 
given by the 2/-direction component of a given vector field. In this example, 
the direction of v needs to not always be identical to that of y and may be 
reverse to that of y. 

Let T X M be the tangential plane at x. Then, we have a loop made of all end 
points of the vectors Y G T X M satisfying 

L{x,Y) = l. (5) 
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This equation is obtained by assuming t = s in Eq. ([3]). Thus, we have a Finsler 
length scale at x such that L(x, Y) = 1. In the case of Eq. (j2J), the condition 
L(x, y) = 1 implies that Jj2i — |v|. Since the length |v| depends on its 

direction, we understand that the Euclidean length yJ2jY? depends on its 
direction on the loop L(x,Y) = 1. To the contrary, the Finsler length of Y 
is constant, which is 1 in the unit of |v|, and hence it is independent of the 
direction on the loop. 

In the example in Eq. (jlj), the Euclidean length of s is given by s\v\, which is 
identical with the Finsler length if | v| = 1, and hence Finsler length s along C 
is implicitly dependent on v along the direction of y. Thus, we find that the 
Finsler length s defined at x on M depends both on the direction of v and on 
the length of v. Since the length of v is dependent on x, the Finsler length 
depends on the position and the direction. 

2.2 Finsler length on triangulated surfaces 

We assume that M is smoothly or piece-wise linearly triangulated such that 
the Euclidean bond lengths are given. On the triangulated surfaces, the ver- 
tices, the bonds, and the triangles are independently labeled by sequential 
numbers. Let N, Nb, Nt be the total number of vertices, the total number of 
bonds, and the total number of triangles. One additional assumption is that 
the bond ij, which is connecting two neighboring vertices i and j, is labeled 
also by velocities or velocity magnitudes and Vji, which are positive. We 
should note that t>y ^ Vji in general and that v ij plays a role of the unit 
of Finsler length from i to j. Thus the bonds are labeled not only by a se- 
ries of integers 1, 2, • • ■ , Nb but also by two series of positive real numbers 
t>i, v 2 , ■ ■ ■ ,vn b and v[, v' 2 , ■ ■ ■ , v' N . It is also possible to assume that v i = v' i . 




Fig. 1. A smooth triangle A in M with a local coordinate (a) at the vertex 1, (b) 
at vertex 2, and (c) at the vertex 3. 

We should comment on a natural assumption that each triangle is labeled by 
a single local coordinate. Indeed, we have three possible local coordinates on a 
triangle, and the total number of coordinates is 3N? for a given triangulation 
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of M. Thus, we chose a set of Nt coordinates from those possible 3 Nt ones. 
The three possible coordinate systems of a triangle A are shown in Figs. QJa)- 
\HyC). In Fig. [TJa), the local coordinate of A is denoted by (xi,x 2 ) such that 
the origin of the coordinate axes coincides with the vertex 1. The velocity 
parameters t>i2 and t>i3 are defined along yi and y2 respectively so that the 
direction of V12 (^13) from 1 to 2 (1 to 3) coincides with the direction of y\ 
(2/2)- Note that the velocity V21 is not included in the coordinate system in 
Fig. [Qa), while it is included in Fig. [Hb). 

Note that the two different velocity Vi and v[ can be obtained from a smooth 
and non-constant vector field v on the surface. Indeed, suppose that v has 
value only on the vertices of the triangulated surface. In such case (see Fig. 
[TJ^a)), if the local coordinate origin is located at the vertex 1, the value V12 
can be obtained from v(l) at the vertex 1. On the contrary, v 2 \ is obtained 
from v(2) in the case of Fig. []Jb). If v is not constant, then v(l) 7^ v(2), and 
therefore we have v 127^21 i n general. 

We firstly define a Finsler function L(x, y) on A in Fig. Ufa) such that 

. yi/«i2 (on the xi axis) , 
L(x,y) = { (6) 

2/2/^13 (on the x 2 axis) . 

Thus the Finsler length L 12 of the bond 12 is given by L 12 = L 1 /f 12 . Indeed, 
L12 = /1 L(x,y)dt = Ji(l/v 12 )yidt = (l/v 12 )yiAt = (l/v^dxv This proves 
L12 = L\/v\2 because dx\ = L\, where L\ is the Euclidean length of the bond 
12. It is also easy to see that L 13 = L 2 /v 13 . 

Secondly, we define a Regge metric [2"T|22f23j on A in Fig. Ufa) such that 

/ 



iJab 



Li 
\F 3 L\ 



F 3 = LiLaCos^g, (7) 



where $3 is the internal angle of the vertex 1. The triangular relation such as 
Li+Lj >Lk is assumed. We should note that Sf=i is not always constrained 
to be 7T. For this reason, a deficit angle tp is defined on the triangle A such 
that ip = J2i=i &i ~ 71 [24J. Equivalently, the internal angle $j is obtained 
from (p such that = $° (1 + (p/rr) , (i — 1, 2, 3), where £f=i $° = it. We 
assume in this paper that the variables $j (i = l, 2, 3), and hence the variables 
Fi (i = 1, 2, 3), are independent to each other for the numerical simplicity. Note 
that g^ b reduces to the Euclidean metric 5 a b if L\ = L-i = \ and -F 3 = 0. 

Note also that the Euclidean edge length Lj is independent of the local coor- 
dinate, while the Regge metric g^ b itself depends on the coordinate. For this 
reason, the discrete Hamiltonian is obtained by using three different g^ b cor- 
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responding to three different coordinates at each triangle in the conventional 
model. 



By replacing L\ by L 12 and L 2 by L 13 in the Regge metric g^ b in Eq. ([7]), we 
have a Finsler metric g^ b such that 



This g^ b gives two different lengths for an edge of the triangle. Indeed, the 
edge length L 12 (= Li/v u ) with respect to the local coordinate in Fig. [T^a) is 
different from the one L 2 \ (= Li/v 2 {) with respect to the local coordinate in Fig. 
QJb) because v\ 2 7^ v 2 \. This is in sharp contrast to the case of the conventional 
Regge metric, where the bond length of A is uniquely given independent of 
the coordinates. 



Fig. 2. (a) A triangle in M with the local coordinate origin is at the vertex 1, and 
(b) a neighboring triangle with the local coordinate origin is at the vertex 2. Two 
triangles share the bond 12. 

Figures [2](a) and E^b) show two neighboring triangles which have a common 
bond 12, where the vertices 1 and 2 are the origins of the two coordinates. In 
this case, both t> 12 and v 2 \ are used to define the model. The Finsler length 
of the bond 12 is given by L 12 in the triangle of Fig. [2]^a), while it is given by 
L 2 i in the other triangle of Fig. [2](b). This is a result of the assumption that 
a triangle should be labeled by a single local coordinate. 

The Finsler area Aa of A is given by the determinant of g^ b such that 



We should note that Aa depends on the local coordinate. However, this does 
not mean that A a in Eq. (Q is ill-defined. In fact, the Finsler length depends 
on the coordinate, it is quite natural that Aa depends on the coordinate. 
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(9) 
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We should note that g^ b can also be obtained from the bi-linear form 



J M 



J2 L 



2 

A- 



L, 



^(Lyvj 2 )yl + (Lyvj 2 )y 2 2 + 2(F 3 /v 12 v l3 ) yi y 2 (10) 



such that 



1 d 2 L\ 



9ab — n 



M 



2 dy a dy b 



This expression implies that g^ b is a (0, 2)-tensor just like an ordinary metric 



9ab{X) 



because y is a (1, 0)-tensor and L 2 M is a function. 



From the Lm in Eq. (TTU|) . we also have the Finsler lengths L 12 and L 13 of the 
bonds 12 and 13 of A in Fig. Ufa). Indeed, we have y 2 = = 0) on the xi 
(x 2 ) axis. Thus, we find from the bi-linear form L\ that the bond length of 
X\ (x 2 ) axis is given by L 12 = Li/v\ 2 (L 13 = L 2 /vi 3 ). We can also start with 
the form Lm in Eq. ( TTDT) . because L\ can be written as L\ = YjabQ^bUayb by 
using g F ab in Eq. 



A surface model with a Finsler metric 




(a) 



(6) 




Fig. 3. (a) A smooth triangle A in M, (b) the image X{A) of Zi by a mapping X 
from M to R 3 , and (c) a unit normal vector no of X(A) and those iij(z = l, 2, 3) of 
the nearest neighbor triangles in R 3 . 

In this subsection, we define a discrete model by introducing the discrete 
Hamiltonian and the partition function. We assume that the surface is em- 
bedded in R 3 by a mapping X : M — > R 3 . A local coordinate is also assumed 
to be fixed on every triangle A in M. 



A discrete Hamiltonian S(X,g) is defined by 



S(X,g)=S 1 + bS 2 , 

S 1 = Ws 1 (A)/A A , S 2 =Ws 2 (A)/A A , 

A * A 



(12) 
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with 

S 1 (A) = (Ll/v 2 12 ) i\ + (L 2 2 /v 2 13 ) t\ - 2 (F 3 /v 12 v 13 ) I ■ £ 2 (13) 

and 



S 2 (A) = (Ll/v 2 12 ) (1 - n ■ m) + (l 2 2 /v 2 13 ) (1 - n • n 2 ) 

- 2 (F 3 /v 12 v 13 ) (n - ni) • (n - n 2 ) . (14) 

The length Lj and F3 = L1L2 cos $3 are the variables of A in M (Fig. [3(a)), 
while the length li and the unit normal vectors n;(i = l, 2) are those of X(A) 
in R 3 (Figs. [3(b), (c)). The coefficients L 2 /v 2 k together with alter the 

surface tension coefficient a(= 1) and the bending rigidity b to be dependent on 
the position and direction on the surface. If is random and hence isotropic, 
the effective surface tension and bending rigidity are expected to be almost 
uniform and not anisotropic. To the contrary, an anisotropic is expected 
to make the coefficients anisotropic. 

The discrete partition function is defined by 

, N 

Z(b) = J Dgj ]JdX t exp[-S(X,g)], (15) 



i=i 



where 



r N B / n b \ 

Dg = J2 ]I dL i ex P ~ X l £ L 2 

C i=l \ i=l / 

JUdFijll dvi exp -\ v V^ 2 • (16) 



i=i " i=i 



/'ni^i^-^i i n ^(^) denotes that the integrations are performed under the 
constraint that the center of mass of the surface X(M) is fixed at the origin of 
R 3 . As mentioned in the previous subsection, we use the variable F in place of 
if. The integration measure JH^dipiexp \ —X ip J2^i \ which is assumed 
in the model of [23], is replaced by /ni=i dFi with a constraint |F 3 | < L 1 L 2 . 
Under this constraint the variable F plays a role of the deficit angle X^=i ^ 
of the triangle. 



The role of the factor exp f — Al J2i=i L 2 j in Eq. ( TIE]) , where Lj is the Euclidean 
bond length, is to suppress the divergence of L. The factor exp (— J2i l/^ 2 ) 
also prevents the variable v from being zero. The distribution of {vi} becomes 
at random and hence defines a random vector field on a surface. If the vector 
field v is given by some other reasons externally or dynamically, this factor 
may be changed. 
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The symbol J2c i n Eq. (TT6j) denotes the sum over all possible coordinates C A 
local coordinate C is fixed on A. In the conventional models such as the model 
in [23], the Hamiltonians S\ and S 2 are defined by including the terms cyclic 
under permutations of three different coordinates of A such that 1 — > 2, 2 — > 3, 
and 3 — )■ 1. A permutation of three different values of C is not a coordinate 
transformation in A. In fact, as C changes from one to another, the discrete 
Hamiltonians S\ and 5*2 change. This is true not only in the conventional model 
but also in the Finsler geometric model. In this sense, C can be viewed as a 
variable just like the triangulation T. However, the dynamical triangulation 
J2r changes the lattice structure. Therefore, J2t cannot be included in J2c- 

We should emphasize that the integrations of the variables F t and v j depend 
on the coordinate C. In each local coordinate of a triangle, only one of the 
two variables, such as vu and v 2 i, is integrated. Analogously in the integra- 
tions /' H^dFi, the variable Fi, which represents one of three different F, is 
integrated. 

2.4 Continuous surface model 

The surface model of Helfrich and Polyakov is defined by a mapping X from 
M to R 3 such that X : M 3 (x 1 ,x 2 ) (->■ X(x 1 ,x 2 ) G R 3 [I]. The symbol 
(xi,x 2 ) denotes a local coordinate of M. 

The Hamiltonian of the model is given by a linear combination of the Gaussian 
potential Si and the extrinsic curvature energy 5* 2 such that 



where b[kT] is the bending rigidity. The matrix g a b{d, b = l, 2) in Si and £2 is 
a Riemannian metric on M , g is the determinant of g a b, and g ab is its inverse. 
The symbol n M in S 2 is a unit normal vector of the surface. We should note 
that S 2 is obtained from Polyakov's action for extrinsic curvature by assuming 
gab = d a X^dbX^ [3]. We here assume that g a b is arbitrary in Eq. ( fT7l) . 

The surface model described by £ in Eq. ([17]) is statistical mechanically defined 
by the partition function 



where S(X, g) denotes that S depends on the variables X and g. The in- 
tegration symbols / Dg and / DX denote the sum over metrics g and the 



S = St + bS 2 




(17) 




(18) 
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mappings X. The model is characterized by the conformal invariance and 
the reparametrization invariance. The first is that the action S remains un- 
changed under a transformation g ab — > g' ab = fg ab for an arbitrary positive 
function f(x). The second is that S remains unchanged under any local co- 
ordinate transformation x — > x'. The transformation x — > x' changes both g 
and X, while the conformal transformation only changes g. 

We simply deform this continuous model by replacing g a b{x) with a Finsler 
metric g a b(x,y), which is a four- variables function. In this new model the 
conformal invariance is apparently preserved even when the factor f(x) is re- 
placed by f(x,y). In contrast, the reparametrization invariance is not always 
preserved, or in other words the reparametrization for x is not always extended 
to the one for x,y(=x). The reason is that the parameter y is not a coordi- 
nate and is only allowed to transform according to a linear transformation 
corresponding to a coordinate transformation of x. However, the Finsler met- 
ric g a b(x, y) is formally a (0, 2)-tensor just like g a b(x) as mentioned just below 
Eq.f TiTj) . Therefore, the continuous actions with a Finsler metric remains a 
scalar and hence are well-defined as action functionals. 

The discrete Hamiltonians in Eqs. (fT2|) - (fl4|) are obtained from the continuous 
actions Si and S 2 in Eq. (TTTj) on the triangulated surface by the replacements 
d x X^^X^-X^, d 2 X^^X^-X^, where Xf denotes the position of the vertex 
i such that t\ = [X^—X^l, £ 2 = |^3 — X± \ (see Figs. Eta)-[3tc)). The derivatives 
in 5*2 in Eq. ffTT|) can also be discretized by din^ — )■ n — n 2 , d^n^ — )■ n — ni, 
where n^(i = 1, 2, 3) are the unit normal vectors shown in Fig. [3]^c). 



3 Simulation results 

3.1 Euclidean model 

The model introduced in Subsection 12.31 is meaningful even in the case where 
the velocity parameter is fixed such that Vij = 1 for all ij and g^ b = 5 a b or 
g^ b = d a X' J 'd b X^. Indeed, the models in those cases are still not always identical 
to the corresponding conventional models because of J2c m Eq. (TTB]) . Therefore, 
in order to see the influence of J2c on ^ ne P nase structure, we study not only 
a non-trivial Finsler geometric model but also the most simple model with 
Vij = l and gf b = S ab . 

The so-called crumpling transition between the smooth spherical phase at high 
bending region and the collapsed phase at low bending region has long been 
studied theoretically and numerically [7|8|9|T0 ,25,26,27j. Both of the phases 
separated by this transition are isotropic in the sense that the surfaces are 
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symmetric under arbitrary three-dimensional rotations. 



To see whether this transition is not influenced by the Finsler geometric treat- 
ment, we firstly study the most simple model, which is defined by 



z ( b ) = e / n dx i ex p [- + bs ?)] > 

C i=l 

s x = (x, - x 3 f , s 2 = (! - * • % 



(19) 



As mentioned above, this model is identical with the conventional model ex- 
cept J2c m Z(b). The Hamiltonians are disctretized on a spherical lattice, 
which is characterized by (N^Nb^Nt) = (10£ 2 + 2, 30£ 2 , 20£ 2 ), where £ is the 
number of partitions of an edge of the icosahedron. 

The canonical Metropolis Monte Carlo technique is used to update the vari- 
ables. The update X' = X+6X is accepted with probability Min[l, exp(— SS')), 
5S' = S"(new) — S'(old), where 5X is a random three-dimensional vector in 
a small sphere. The radius of this sphere is fixed to a constant to make the 
acceptance rate of X approximately 50%. 

The sum over coordinates J2c m ls performed as follows: On a triangle 
A, the coordinate is characterized by its origin, and hence A has only three 
possible coordinates. Therefore, the current coordinate of A is randomly up- 
dated to one of the two remaining coordinates. In this update of C, Si and S 2 
change from the expressions in Eqs. f )13p and fll4p to those cyclic under the 
permutations 1 — » 2, 2 — > 3, and 3 — > 1, and so on. 




Fig. 4. Three configurations of coordinate axes in two neighboring triangles. Neither 
the bending energy nor the bond potential is given to the bond 12 in the coordinate 
of (a), while both of the energies are redundantly given to the bond 12 in the 
coordinate of (c). 

No constraint is imposed on the update of C. As a consequence, the configu- 
ration like the one in Fig. 11(a) appears, where the bond 12 shares neither S\ 
nor S^. Thus, the random update of C makes the configuration non-uniform in 
the sense that the distributions of S\ and S 2 are non-uniform. Therefore, it is 
non-trivial whether the phase structure of the conventional model is influenced 
by J2c in Z. 
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1.29 



1.31 b 1.33 





-> — i — 1 — i — 1 — i — 1 — r 

5^_N=15212 




N= 12252 




N=8412 
_jV N=4842 






Uvf N=2562 


1.29 


1.31 b 1.33 




Fig. 5. (a) The mean square size X 2 vs. b, (b) the variance Cx^ vs. 6, and (c) a 
log- log plot of the peak C^f x vs. N. The solid lines connecting the data in (a),(b) 
are drawn to guide the eyes. The straight line in (c) is drawn by fitting the data to 
Eq. ([22]). 



The mean square size X 2 is defined by 



X 



X 



X 



N 



(20) 



where X is the center of mass of the surface. The large errors in X 2 reflects 
a large fluctuation of X 2 as a result of the crumpling transition between the 
smooth phase and the crumpled phase (Fig. 0(a)). The variance C x ^ of X 2 
defined by 

^ = ^((X 2 -(X 2 )) 2 ) (21) 

can reflect the phase transition. The peak at 1.303 in Fig. Efb) indicates 
the existence of the transition. The peak values C™2 X scales against iV; the 
straight line in Fig. E(c) is drawn by fitting the data to 



/imax 



N a , 



a 



1.18 ±0.07. 



(22) 



The obtained exponent a 
order. 



1.18(7) indicates that the transition is of first- 



We find from the plot of the series {X 2 } in Fig. E](a) that the smooth and 
crumpled phases are clearly separated. This series is obtained at the transition 
point b = 1.304 on the N = 15212 surface. Plots of {X 2 } similar to the one 
in Fig. [6](a) are obtained on the N = 12252 and N = 8412 surfaces, though 
they are not depicted. In order to have the Hausdorff dimension H defined by 
X 2 ~ N 2 I H , we calculate the mean values of X 2 in the smooth and crumpled 
phases independently from the series {X 2 }. On the surfaces N > 8412, we 
use a series {X 2 } at the transition point like the one in Fig. [6](a), while on 
the surfaces iV<4842 we use two different {X 2 } obtained in the smooth and 
crumpled phases. Figure MJo) shows the results X 2 vs. N in a log-log scale. 



12 




Fig. 6. (a) The variation of X 2 against MCS at the transition on the N = 15212 
surface, (b) log-log plots of X 2 vs. N in the smooth and crumpled phases at the 
transition point, (c) S\/N vs. b. 

From the slope of the fitted lines, we have 

H sm = 2.20 ± 0.28 (smooth), H cr = 2.57 ± 0.95 (crumpled). (23) 

These values are comparable to the results H sm = 2.02(14) and H cr = 2.59(57) 
of the conventional model within the errors |28j . 

We see the expected relation S±/N— 1.5 in Fig. 0(c). This implies that the 
equilibrium configurations are correctly obtained under the treatment J2c m 
Z. 




Fig. 7. (a) The bending energy S^/iVy vs. b, (b) the specific heat Cs 2 vs. b, and (c) 
a log-log plot of the peak values C<? ax vs. N. The fitting is performed by using the 
largest three data in (c). 

The bending energy S%/Nt, the specific heat 

^ 2 = |((S 2 -(S 2 >) 2 ), (24) 

and the peak values C™ x are plotted in Figs. C3a)-[7tc). By fitting the largest 
three values of C^ ax , we have 

C£ ax ~ N a , o — 1.02 ± 0.16. (25) 
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The exponent a — 1.02(16) also confirms that transition is of first-order, and 
the value is in good agreement with a = 0.93(13) of the conventional model 



3.2 A Finsler geometric model 

In this subsection, we study a non-trivial model, which is defined by 

N T N 



Z(b) j n>l j\\d,>, j J! dXi exp [- (Si + bS 2 + X 3 S 3 + A 3 S 4 )] , 

i=l i=l 
N B I N B \ N B / 

/ d 9=t,[ n^ ex p -^x>. ? / n^ ex p - a «E!a 

C i=l \ i=l / i=l \ i 

5i and S 2 = Eqs. (HI,®, flU, (26) 

N T 

s 3 = ^piAi, s 4 = \Pi - PjV 

i=l ij 

In this model, the variable (p for the deficit angle is suppressed for simplicity, 
and hence the variable F in Eq. f|T6|) is only given by F 3 = (L\ + L\ — L\) /2. 
In this case the metric g^ b in Eq.(J7J) is identical with the conventional Regge 
metric [2T1I221I23] 

The symbol pi in S 3 is a scalar field on M, which is a conjugate variable to 
the surface area A i: and S4 is the interaction term. m 5*4 is the sum over 
all nearest neighbor triangles i and j. The coefficients A 3 and A 4 are fixed to 
A 3 = A 4 = 1. The variable pi and the terms S3 and S4 are not always necessary 
but they are introduced to take an in-plane deformation into account. In this 
model, p interacts with the surface through the coupling piAi. To the contrary, 
a constant scalar field p has no explicit interaction with the surface. If A$ is 
constant, then the field p and hence both 5*3 and 5*4 are independent of the 
surface geometry. 

The effective Hamiltonian S' including the measure terms is given by 5" = 
Si+bS2+J2iL^+J2i^i+S 3 +S4, where Lj is the Euclidean bond length. In this 
expression, we replace 1/fj with the inverse velocity Vi = 1/vi for numerical 
simplicity. The variables summed over in the partition function are X, L, v, 
and C. One MCS consists of iV updates of X, Nb updates of L, Nb updates 
of Vi, and Nt updates of C. 

The variables L is updated such that V = L+5L(> 0) with random numbers 
SL G [—0.5, 0.5]. In this update, V is constrained to satisfy the triangle equal- 
ities. The inverse velocity Vi(i = 1,2) is updated such that v[ = Vi + 5vi(> 0) 
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with a random number 8vj G [—0.5, 0.5]. None of the variables Vi is updated on 
the bond 12 of the configuration in Fig. BJa), one of is updated in Fig. BJb), 
and both of the variables are updated in Fig. BJc). The constraint Vi < 1 is 
imposed. Without this constraint, the acceptance rate for the update of local 
coordinate C remains very small (10% ~ 15%), while it remains 40% ~ 50% 
under the constraint. 




Fig. 8. (a) The mean square size X 2 vs. b, (b) the variance Cx^ vs. b, and (c) a 
log-log plot of the peak values C^f* vs. N, where A3 = A4 = 1. 

The mean square size X 2 against b in Fig. [Ha) is almost identical to that 
of the Euclidean model in the previous subsection. The variance Cx? and the 
peak value C^f* shown in Figs.[8](b) and[8]^c) are also almost identical to those 
of the previous subsection. The straight line in Fig. [8]^c) is drawn by fitting 
the data to Eq. (12"2"I) . and we have a = 1.28 ± 0.03. This value is identical to 
the one in Eq. (1221) within the error. 




Fig. 9. (a) The bending energy S^/A^ vs. b, (b) the specific heat Cs 2 vs. b, and (c) 
a log- log plot of the peak values C<? ax vs. N, where A3 = A4 = 1. 

The bending energy S2/NB, the specific heat Cs 2 are almost identical with 
those of the Euclidean model (Figs. [^a),(b)). The scaling of C!? ax predicted 
in Eq. ( 12"5"|) gives the exponent a = 0.95(6), which is almost comparable to 
that of the Euclidean model. The Hausdorff dimension H is calculated from 
the series of {X 2 } at the transition, and the results are H sm = 2.12(20) and 
H cr = 3.0(10). The result H sm = 2. 12(20) is comparable the one in the previous 
subsection, while H cr = 3.0(10) is slightly larger the corresponding result in 
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the previous subsection. However, we see no difference in the phase structures 
between the Finsler and conventional models. This implies that the Finsler 
geometric treatments including J2c are well-defined. 

We performed the simulations for a model with the variable F, which corre- 
sponds to the deficit angle (p of the triangles in M. This model is identical with 
the one in this subsection except the variable F. The results are consistent 
with those of the conventional model just like the models shown in this and 
the preceding subsections. 



3.3 Anisotropic surface model 

In this subsection, we see that a tubular surface is obtained under a constant 
vector field on the surface. Let g^ b in Eq. ([7]) be g^ b = 5 a b and \v x \ = t>i 2 = v, 
\vy\ — v i3 — 1) where (x,y) is a local coordinate along v on M, then we 
have 5*2 = (1/2) / dxdy [(l/v)(d x n^) 2 + v {d y n^) 2 ]. Thus we have an anisotropic 
bending rigidity such that b x = b/v and b y = bv. If v<l, then we have b x >b y , 
and consequently the surface becomes smooth (wrinkled) in the x direction (y 
direction) in a certain range of b. In this model, a tubular surface is expected 
at sufficiently large or small v, although none of the parameters b x , b y can be 
exactly 0. 

Thus we have anisotropic bending rigidities b x and b y if the vector field v is 
constant. Since a vector field v on M corresponds to the one on the surface 
X(M) C R 3 , we assume a constant in-plane tilt order a at the center of 
each triangle. This in-plane variable a becomes a vector field on X(M) and 
corresponds to v in M. 

The variable a is of unit length \a\ = 1 and has a value in R 3 . The direction 




(a) (6) 



Fig. 10. (a) The x-direction component c?j defined by Eq. (|27p . a unit normal vector 
nj, and the canonical basis vector e x = (1,0,0), of the triangle i, (b) the variables 
<7i and $j of the triangles i and j, and a unit tangential vector of the bond ij. 

of &i of the triangle i is defined by the projection of e x = (1, 0, 0) G R 3 on the 
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triangle plane (Fig. HPT a)) such that 



0i = i 7 \ — r> I 27 ) 

\e x - (e x . • Hi) ni| 



where is a unit normal vector. As the surface shape varies, not only nj but 
also <7j varies. We use the word " constant" in the sense that <7j is defined by a 
constant vector e x = (l, 0, 0). This a, represents a constant in-plane tilt order. 

The variable &i plays a role of v on the surface, and therefore the component 
&i ■ tij is used to define the bending rigidity at the bond ij, where is a unit 
tangential vector along the bond ij. The effective bending rigidity b{j is given 
by bi2 = bvi 3 Vi2 and 613 — bv nVii, where v±2 and t>i3 are the vectors along the 
bonds 12 and 13 in Fig. Ufa). However, this b^ becomes singular if Vij — 0. For 
this reason, we simply define i>y by multiplying the integer N v so that Vij has 
an integer value in {1, 2, • • ■ , N v }: 



v 



ij 



l + [(Tij], (Tij = N v \B { ■ tij\ , (easel), (28) 



where [x] represents the integer Max{n e Z\n < x}. Consequently, the mini- 
mum (maximum) value of the effective bending rigidity becomes b (N v b). 

Another possible discretization of 6y is given by modifying 6^ from the naive 
discretization bvi = bv^v i 2 l and 613 = bv nv^i to b\2 = bvi 2 l and 613 = fru^ 1 , and 
by defining v^ such that 



-' ) 1 (l*i'td<0.5) 

= < , (case 2), (29) 

&o >0.5) 



where b (< 1) is a constant to be fixed. It is possible to define v~^ in Eq. (j2"9~j) 
by using (1/2) (oj + c?-,) ■ t^. In this case, % is uniquely defined on the bond 
ij (vij = Vji), and hence it is possible to assume the conventional discretization 
for 5*2, where J2c * s not included in the partition function Z. However, we use 



v 



1.1 



1 in Eq. (!29l as a demonstration, where J2c ^ s included in Z. 



The partition function Z, the Gaussian bond potential Si and the bending 
energy S 2 are given by 
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z ( b ) = e / n dx i ex p t- (Si + hs ^ > 

C i=l 

^^(Zi), S 2 = £S 2 (Z\), (30) 
5 1 (Zi) = (X 2 -X 1 ) 2 + (X 3 -X 1 ) 2 , 

S^Zi) = — (1 - n • n x ) + — (1 - n ■ n 2 ) , (case 1), 
S 2 (A) = (1 - n ■ ni) + v[~3 (1 - n ■ n 2 ) , (case 2), 

where the surface tension coefficients a\% = v i 2 v j^ 1 and ai3 = v i$v i 2 of Si(A) are 
fixed to 1 for simplicity. The unit normal vectors no, ni, n 2 in ^(Zi) are those 
shown in Fig. [3](c). Due to the coefficient %(= fisf^ 1 , ■ • • , v {3 1 ) in 82(A), we 
have the effective bending rigidity bbij as mentioned above. 



(a) (b) (c) 




(d) (e) (f) 

Fig. 11. Snapshot of surfaces of size iV = 10242 of case 1 model obtained at (a) 6=1, 
(b) 6 = 2, and (c) 6 = 5, where ^ = 10. Snapshots of case 2 model obtained at (d) 
(6 , 6) = (0.2, 3.5), (e) (6 , 6) = (0.2, 7.5) and (f) (6 Q; 6) = (0.2, 15). 

The integer N v is fixed to N v = 10 in case 1, and the parameter 60 is fixed to 
6o = 0.2 in case 2, while b is varied in the simulations. Snapshots in Figs JTlT a)- 
(f) show that tubular surfaces are obtained in both case 1 ((a)-(c)) and case 
2 ((d)-(f)). Since a vector field on a sphere has singular points, the variable cij 
becomes singular on the surface. Under the condition of Eq. (1271) . we expect 
that there appear two singular points. Indeed, these two points can easily be 
seen at two terminal points of the surface in the snapshots at relatively small 
b in both case 1 and case 2. 
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4 Summary and Conclusion 

We have numerically studied a Finsler geometric (FG) surface model, which 
is defined by a mapping X from a two-dimensional Finsler space M to R 3 . 
The FG surface model is obtained from the model of Hefrich and Polyakov 
(HP) for strings and membranes by replacing the Riemannian metric with a 
Finsler metric g^ b (x,y). In this sense, this model is an extension of the HP 
model. By discretizing this continuous FG model, we have a discrete FG model. 
The discrete Finsler length depends on the direction of local coordinates on 
the surface, where a new fictitious variable C is introduced. The variable C 
represents a coordinate on the triangles. 

To check whether the FG model is well-defined, we firstly study the most 
simple model without the variable v, where gab( x ) = $ab- in this model, the 
crumpling transition between the crumpled and smooth phases is found to 
be of first-order and is identical to that of the conventional model. It is also 
found that a FG surface model with random vector field v has the same phase 
structure as the conventional model. To demonstrate that a tubular surface 
appears, we define g^,{x) by using a constant vector field v. The Monte Carlo 
results show that a tubular surface is obtained in a certain range of the bending 
rigidity b. Thus, we confirm that an interaction between the surface and an 
anisotropic vector field is a possible origin of the surface anisotropy. 

It is interesting to study an anisotropic surface model with a vector field v 
which is not constant or random but dynamically changed. A variable a can 
be defined at the vertices of triangles as a tilt such that its in-plane component 
plays a role of a vector filed on X(M), and interactions between the tilts a 
can also be assumed. The fictitious variable C is summed over in the partition 
function of the models in this paper, however, it can be summed over in 
the discrete Hamiltonians as in the conventional treatment of surface models. 
Studies on the FG surface model should be performed more extensively. 
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